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Abstract 

We conjecture chronology is protected in string theory due to the condensation of 
light winding strings near closed null curves. This condensation triggers a Hagedorn 
phase transition, whose end-point target space geometry should be chronological. 
Contrary to conventional arguments, chronology is protected by an infrared effect. 
We support this conjecture by studying strings in the O-plane orbifold, where we 
show that some winding string states are unstable and condense in the non-causal 
region of spacetime. The one-loop string partition function has infrared divergences 
associated to the condensation of these states. 
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1 Introduction 



The possibility of time travel has been one of the theoretical physicists' favorite puzzles. 
Of particular interest is time traveling to the past, since it raises many causality paradoxes. 
Although our experience tells us such voyages are unlikely to be achievable, there is no 
rigorous proof of their impossibility. In an attempt to address this problem, Hawking put 
forward a chronology protection conjecture [1], extending earlier results by Tipler [2]. His 
argument was based on the behavior of quantum field theory in the presence of closed 
null curves, since in a spacetime with a Cauchy surface any future causality problems 
must be preceded by the formation of one such curve. Hawking argued that the one-loop 
energy-momentum tensor becomes very large near a closed null curve, therefore causing 
a significant backreaction in the geometry. The outcome would be a naked singularity or 
a gravitational evolution that prevents formation of the closed null curve itself. Either 
way, the mechanism that protects chronology is an ultraviolet effect. Since the ultraviolet 
behavior of string theory is different from that of quantum field theories, it is questionable 
that the above mechanism will protect chronology. 

Theories with extended objects provide privileged probes of spacetimes with closed 
causal curves since the extended objects can wrap around them. As mentioned above, 
closed null curves are particularly important because causality problems first arise when 
such curves are created. In the framework of string theory, there are winding string states 
that can become light just before they wrap a closed null curve, when their proper length 
is of the string scale. A phase transition can then occur, analogous to the Hagedorn phase 
transition where winding string states become massless for a compactification circle with 
radius of the string scale [3, 4]. It is therefore plausible that light winding states play a 
role in protecting chronology. In this work we establish this analogy in detail by studying 
a toy model, where we are able to show that such phase transition does in fact occur and 
it is associated to the condensation of string fields in the causally problematic region. We 
expect this mechanism to hold in general, which leads us to the following conjecture: closed 
null curves do not form in string theory because light winding states condense, causing a 
phase transition whose end-point target space geometry is chronological. This mechanism 
for protecting chronology is truly stringy and is an infrared effect, in sharp contrast with 
Hawking's proposal. 

The toy model we shall consider is an orbifold of three-dimensional Minkowski space 
introduced in [5], called the O-plane orbifold. The orbifold quotient space is a good labo- 
ratory to study the issue of chronology protection because it is free of singularities and it 
develops closed causal curves. There is, however, a truncation of spacetime that is free of 
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closed causal curves. It is the main point of this work to show that some winding string 
states, which condense in the pathological region of spacetime, give a possible mechanism 
for the protection of chronology. For other works of strings on orbifolds with closed causal 
curves and on closed string backgrounds with bad chronology see [6-47]. 

Since the orbifold closed causal curves are topological, it is questionable whether the 
conjectured mechanism is generic, in the sense that it might not work for non-topological 
closed causal curves dictated by the Einstein equations [48]. However, when fermions obey 
anti-periodic boundary conditions around the closed curves, we believe the mechanism is 
general because the appearance of unstable winding states relies only on the existence of 
closed curves with very small proper length. 

We start our discussion in section two by giving a detailed review of the O-plane 
orbifold. In particular, we analyze its causal structure, wave functions and quantum field 
theory divergences. We find infrared divergences associated to the condensation of some 
particle states and show that the theory is ultraviolet divergent because of the presence 
of the non-causal region. In section three we consider strings in the O-plane orbifold. 
Again we analyze the string wave functions and compute the one-loop partition function, 
which when appropriately interpreted has only infrared divergences. These divergences 
are associated to the presence of unstable string states, whose condensation is conjectured 
to protect chronology in section four. We give some concluding remarks in section five. 
Appendix A contains some technical results regarding wave functions and corresponding 
Hilbert space measure. In appendix B we give an alternative derivation of the particle and 
string partition functions using the path integral formalism. 



2 Review of O-plane orbifold 

In this section we shall review some basic facts about the O-plane orbifold introduced in 
[5]. This is an orbifold of three-dimensional Minkowski spacetime. To define it consider 
the light-cone coordinates ( X j X j X + ), normalized such that the flat metric is given by 
ds 2 = —2dx + dx~ + dx 2 . The orbifold group element VI = e K is generated by the Killing 
vector 

k = 2ni (RP_ + A J) , 

where 

d d d 
tj — x+ ~~r — %~r\ — T~ 1 ^ 



dx dx + ' dx~ ' 

are, respectively, the generators of a null boost and a null translation and R and A are 
constants. The orbits of this Killing vector field can be found in figure 1. The orbifold 
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Figure 1: Orbits of the Killing vector field k projected in the XX~-plane of the covering space. 
The surface k 2 = divides spacetime in a good region and a bad region. The dashed straight 
lines are timelike geodesies connecting image points, projected in the XX~-plane. 
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Since the generator of the O-plane orbifold involves the time direction one might wonder 
about the causal structure of quotient space, and in particular about the existence of closed 
causal curves. For this purpose it is convenient to introduce a new set of coordinates 
y = (y~ ,y,y + ), in which the orbifold action becomes trivial. These are defined by 

x~ = y~ 

x = y-^(y~) 2 

E 2 

x + = y + -Eyy +-g-(j/ ) 3 , 
where E = A/R. Now the orbifold identification becomes simply 

y~ ~ y~ + 2-kR 

and the metric 

ds 2 = -2dy-dy + + 2Ey(dy~) 2 + dy 2 . 
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The surface y = divides spacetime in two regions, one with k spacelike (y > 0) - the 
good region - and the other with k timelike (y < 0) - the bad region. The region with 
y < is expected to be pathological. Indeed, although there are closed timelike curves 
passing through every point in this space, these curves always enter the y < region as 
shown in figure 1. Moreover, in the covering space there are points which are light-like 
separated from their image. These points lie on the so-called polarization surfaces [49]. 
Null curves joining a point and its image become closed null curves in quotient space and 
are dangerous in quantum field theory because they are responsible for extra divergences, 
as we shall review in subsection 2.3. 

2.1 Particle dynamics 

To understand dynamics in the O-plane orbifold, we start by analyzing the motion of point 
particles. Clearly, particle trajectories are straight lines in covering space. However, we are 
interested in studying the problem from the view point of an observer using ^-coordinates, 
which reduces to considering the following Lagrangian 

£ = -2y+y- + 2Ey{y-f + y 2 ■ 

Dots are derivatives with respect to an afline parameter. We now introduce the canonical 
conjugate momenta (p_,p,p + ), of which the first and third components are conserved. This 
is because the quotient space preserves only two of the six isometries of three-dimensional 
Minkowski space, namely those generated by d + = d/dy + and cL = d/dy" . The problem 
then simplifies to the motion of the particle along direction y, with an effective Hamiltonian 
given by 

H = p 2 + V(y) , V{y) = -2Ep 2 + y-2p + p_ . (1) 

The motion is thus determined by a linear potential. The constant force associated to this 
potential is clearly an inertial force, since the ^/-coordinates describe an accelerating frame. 

Now comes an important point. The Hamiltonian, which is a constant of motion, is 
the momentum squared of the particle 

n = v ^ = -M\ (2) 

and is therefore fixed by the particle's mass. The classical turning point yo for a particle 
of light-cone energy p + and Kaluza-Klein momentum p_ is then found by setting p = in 
the previous equation 

M 2 = 2p + p^ + 2Ey pl . (3) 
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Figure 2: Effective potential V(y) and conserved energy — M 2 for particle motion along direction 
y. Shown are the cases of M 2 > (left) and M 2 < (right) for p_ = 0. The turning point yo is 



Consider first a particle of M 2 > and p- = 0. It is clear from figure 2 that such 
a particle will never enter the bad region y < 0. Since the slope of the potential gets 
steeper as p + increases, the higher the light-cone energy of the particle, the closer it gets 
to y — 0. In the low energy limit p + — > 0, the particle is pushed away to y — +oo. On the 
other hand, particles of M 2 < and p_ = always cross into the bad region and, for low 
energies, reach very far out into this region as can be seen from figure 2. 

For particles with p_ ^ the potential is displaced along y, so it no longer goes through 
y — 0. However, the above qualitative behavior is the same: particles with M 2 > are 
repelled away from the bad region and for low energies are pushed away to y — +oo; 
particles with M 2 < are also repelled but for low energies reach far out into the bad 
region. 

2.2 Single particle wave functions 

To quantize the particle, we follow the usual prescription of replacing p^ — > — iD^ in the 
mass-shell condition (2). The single particle wave functions (p(y~ ,y,y + ) obey the Klein- 
Gordon equation 



To find the solutions to this differential equation we first consider a complete basis of 
functions in the orbifold. The details are presented in the appendix A.l; here we state the 
necessary results. In addition to the orbifold periodicity < y~ < 2nR, we consider a 



defined by V(yo) = —M 2 . Particle states with p 2 < are allowed quantum mechanically. 
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Figure 3: Profile of a wave function in the y-direction. 



large box defined by < y + < L + and — L/2 < y < L/2. The wave functions 

\K{p + )\ 1 '* 



y ,m\ 



27rRL + p(y ,p + ) 



Ai(z) e 



(4) 



form a complete normalized basis and are eigenfunctions of the Laplacian in the O-plane 
orbifold. The Kaluza-Klein momentum is quantized as p- = m/R and the other quantum 
numbers are the light-cone energy p + and the classical turning point y$. The function 
Ai(z) is an Airy function with argument given by z 3 = K(y Q — y) 3 and 



K = K{p + ) = 2Ep 



+ ' 



p(yo,p+) = - 

71 



K(p+) ( sgn(K) --Vo 



1/2 



(5) 



A y-plot of a wave function is given in figure 3. To sum over states in the Hilbert space 
the appropriate integration measure is given by 



E 



2n 



dp+ / \ dy p(y ,p + ) . 



On-shell wave functions are eigenfunctions of the Laplacian with eigenvalue M 2 and 
therefore their quantum numbers satisfy the mass-shell condition (3). To define the evolu- 
tion of a given field configuration, we shall consider an initial value problem defined on a 
surface of constant y + . The normal h to this surface has norm given by n 2 = —2Ey. Thus, 
in the good region y > this surface is space-like and the initial value problem is standard. 
However, for y < we are specifying the "initial" data on a time-like surface. We shall 
interpret this "initial" data as a boundary condition. This is satisfactory to understand 
the evolution of perturbations in the good region. In appendix A. 2 we construct a basis 
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Figure 4: Dispersion relation for particles with M 2 > and m = 0. 

of functions on a surface of constant y + . The light-cone evolution of each mode will be 
determined by its light-cone energy fixed by the on-shell relation. In the appendix we 
show that, in order to be normalizable on a surface of constant y + , the modes must have 
K G R. In particular, for K < 0, we have unstable modes with p + purely imaginary and 
yo complex. 

To understand the subtleties of matter propagation in the O-plane orbifold, we consider 
the cases of massive particles (M 2 > 0) and tachyons (M 2 < 0) separately, both for m — 0. 
Massless particles can be included in the case of massive particles by introducing an infrared 
regulator. 

For massive particles, the behavior of the wave function can be easily understood from 
the corresponding classical effective potential in (1). The turning point y is now the point 
where the wave function turns from oscillatory to evanescent (see figure 2). For on-shell 
particles, the mass formula (3), with p_ = 0, defines the light-cone energy for a given y . 
This dispersion relation is plotted in figure 4. As expected from the classical behavior, in 
the low energy limit we have y — > +oo and the wave function vanishes everywhere. In 
appendix A. 2 we show that the set of on-shell wave functions with y G R is a complete 
basis of normalizable functions on a surface of y + constant. Thus, from figure 4, we see 
that to include modes with y < we must allow for imaginary p + . Then we are effectively 
dealing with the classical potential V(y) with p 2 < and the wave functions oscillate in 
the bad region. Since the light-cone energy is imaginary, these modes can grow in time. A 
generic field configuration will then grow in light-cone time y + , mainly in the bad region. 
We conclude that massive particles have unstable modes in this spacetime. If the bad 
region were to be excised, we would expect to be left only with well-behaved modes. 

Next consider a tachyon, also with m — 0. As shown in figure 5, for p + real the turning 
point lies at yo < and the wave function oscillates according to the classical trajectory in 
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Figure 5: Dispersion relation for particles with M 2 < and m = 0. 



figure 2. To allow for states with yo > one must consider imaginary p + . These unstable 
modes oscillate in space for —oo<y<yo and therefore reach far out into the good region. 
Contrary to massive particles, in the low energy limit the tachyon wave function is a plane 
wave throughout space, which is the usual result for tachyons in Minkowski space. 

Particles with non-vanishing Kaluza-Klein momentum have a similar qualitative be- 
havior. In particular, the condition for having on-shell stable states with real yo can be 
obtained from the mass formula (3) 



where R 2 (yo) = 2Ey R 2 is the proper radius squared of the orbifold compact direction, 
computed at the turning point yo. When this condition is not satisfied, p + becomes com- 
plex and the associated wave functions are not normalizable on a surface of constant y + . 
However, in appendix A. 2, we show that there are other unstable on-shell states, with p + 
purely imaginary and yo complex, which are normalizable. 

2.3 Quantum field theory and ultraviolet catastrophe 

To compare with the string theory results presented later on in this paper, we study the 
usual quantum field theory pathologies in the O-plane orbifold due to the presence of closed 
causal curves. These are the ultraviolet problems that led Hawking to conjecture that a 
strong gravitational backreaction would protect chronology. 

Consider first the existence of polarization surfaces, arising from light-like separated 
images. Since the covering space is Minkowski space, it is simple to compute the relativistic 




(6) 
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interval between image points 

Al = \\x-xJ 2 = 8E{nnR) 2 (y-^{nnR) 2 ^) . 
In particular, points located at 

V = Vn = ^ (nnR) 2 , 

are null separated from their n-th image. It is important to note that the corresponding 
closed null curve intersects the y = surface, so it is possible that cutting off the y < 
region cures any pathology arising from the polarization surfaces, as proposed in [13, 18]. 

For a scalar field of mass M the two-point function can be written in Euclidean 
Schwinger parameterization as 

f ~MA w r+OQ A] \ / / A 2 \ 

<^>^'» = ? W = ? I iww^A-h M2 - it) • 

where = \\x — x^\\ 2 and the trajectory winding number w runs from — oo to +oo. 
When x lies on the n-th polarization surface, the two-point function has extra divergences 
as x' — > x, coming from the w = ±n terms in the sum. These are ultraviolet divergences 
that are not renormalizable [49]. Likewise, the expectation value of the energy-momentum 
tensor diverges at the polarization surfaces. Again, these problems could be cured if one 
excises the bad region. It is the main point of this paper to suggest that a truly stringy 
phenomena in fact excises such region, hence protecting chronology. 

Next consider the one-loop contribution to the vacuum energy, written as 



^ -Tr exp - t -(M 2 + Pfl p^ 



where the trace is over off-shell single particle states. This is written in first quantized 
language so as to compare it later with the string torus amplitude. Choosing the basis 
\p + ,yo,m) to perform the trace and Wick rotating / — > —il, p + — > ip + , we arrive at the 
expression 



Z = t J 21 22 ^ J dp+ dyo p ( yo > W+) exp _ 2 I ~ P+ R + 2EyoP + 

One could question the validity of this Wick rotation since the final expression clearly 
has large I divergences that were absent in the original one. The alternative would be to 
choose a state dependent Wick rotation depending on the sign of (M 2 +p^ l p fl ). However, 
in string theory the requirement of modular invariance forces us to treat all modes equally 
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[14]. The chosen Wick rotation is physically reasonable because massive particle states 
that propagate in the good region, i.e. states with y > 0, give a finite contribution to the 
partition function. On the other hand, states with y < will originate infrared divergences 
provided p + is large enough, for both massive and tachyonic particles. 

Restricting to states with y > and for large /, the p + integral can be done using the 
saddle point approximation. The / integral is then dominated by an exponential with the 
following argument 



M 2 + 



m 



2 



R 2 (yo)J ' 

Hence, the condition for this integral to be infrared convergent is equivalent to the condition 
for the existence of on-shell stable states with y positive, as can be seen from (6). For 
massive particles, this condition holds for all on-shell states that oscillate in the good 
region, which have yo > 0. On the other hand, for tachyonic particles, the above condition 
is not always satisfied and therefore the partition function becomes infrared divergent. This 
is related to the existence of on-shell unstable states. These states oscillate in both good 
and bad regions and therefore their condensation would change the whole of spacetime. 

To analyze the ultraviolet behavior of the partition function it is useful to consider the 
expression obtained from the path integral formalism. We report on this computation in 
appendix B and on its derivation starting from the above canonical result in appendix A.l. 
Here we simply give the final result 

dl 1 / I ^ 2 (2nwR(y)f\ 



/r°° dl ' 



(2tt/) 3 / 2 1 \ 2 21 



where y is now the particle's average position along this direction and w is the winding 
number of closed trajectories. The w = term gives the usual renormalizable ultraviolet 
divergence from the region I — > 0. The other terms are ultraviolet finite for y > 0, diverge 
for y — > and acquire, by analytic continuation, an imaginary part for y < 0. These are 
extra ultraviolet divergences that are absent in string theory. 



3 Strings in the O-plane orbifold 

We now move on from particles to strings. We work with the bosonic string in critical 
dimension, adding the transverse space M 23 to the O-plane orbifold. Since these directions 
are mere spectators we include their contribution explicitly only when necessary. We work 
in units where a' — 1. 
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When considering strings the novelty is the appearance of twisted sectors, which cor- 
respond to strings winding around the ^"-direction. As we shall see, these winding states 
play a crucial role in protecting chronology. 

Start with the string Lagrangian 



-Y+Y~ + EY{Y~f + — + (Y + Y(Y-y - EY(Y~)' 2 - — 

2 2 



(7) 



where dots and primes denote derivatives of the worldsheet fields (Y~ , Y, Y + ) with respect 
to the worldsheet coordinates r and a, respectively. The coordinate a runs from to 
ii. Consider the center of mass motion of a closed string that winds w times around the 
compact direction. The embedding functions have the form 

Y-(r,a) = 2Rwa + y'(r) , Y(r,a) = y(r) , F + (r, a) = y + (r) , (8) 

from which we derive the string center of mass conjugate momenta 

p- = -\ {y + - iEyy-) , p = \y> p+ = ~\y~ ■ ( 9 ) 

Again, only p + and p_ are conserved. As for the particle case, the motion along y is 
determined by the effective Hamiltonian 

U = P 2 + V(y) , V(y) = -2E (p 2 + - (wR) 2 ) y - 2p+p_ . (10) 

We see that there is an additional linear contribution to the potential from the winding of 
the string. Its interpretation is simple. Since the proper length of the compact direction 
increases with y, it costs energy for a winding string to increase its value of y. Defining 
the conserved quantity A = — Tl, the string center of mass turning point yo, where p — 0, 
satisfies 

X = 2 P+ p_+2Ey (pl-(wR) 2 ) . (11) 

For the particular solution (8) the Hamiltonian in (10) generates translations along the 
worldsheet time and it is related to the Virasoro zero-modes by Ti = 2(L + L ). This is 
zero by the classical Virasoro constraints. As we shall see in subsection 3.2, in the quantum 
theory 2(Lq + Lo) = Ti + 2(iV + N) = 4, where N and iV are the usual number operators. 
For each string state it is therefore necessary to consider both negative and positive values 
of A = 2(n + n — 2). The corresponding center of mass classical dynamics reduces to the 
motion of a particle of constant energy in the potential in (10). We analyze it below, 
together with the associated wave functions. 
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3.1 Winding states dynamics and wave functions 

We have just seen that on-shell classical string states satisfy the relation A = —Ti, which 
using (10) yields 

-p^ = A + 2Ey {wRf . (12) 
This constraint defines the wave equation for the string center of mass 

( - 2«9 + «9_ - 2Eydl + d 2 )<P = M 2 (y) , 

where 

M 2 {y) = \ + 2Ey(wR) 2 . (13) 
This Klein-Gordon equation can be derived from the Lagrangian 

C = d<f)d0 + M 2 (y)(f)0 

for a complex scalar field associated with strings of winding number ±w. The major 
difference to the particle case is that this string mass depends on the ^-coordinate and 
therefore it is not conserved. From a pure field theoretical view point one would then 
expect that something special happens for M(y) = 0. We shall come back to this point in 
section 4. 

Solutions to the above differential equation have the same form as for the particle case. 
The only difference is that K in the expression (4) for the basis of functions and in the 
corresponding definitions (5) is now given by 

K(p + ) = 2E (pi - (wR) 2 ) . 

On-shell states satisfy the string dispersion relation (11). As for the particle case, on-shell 
wave functions, which are normalizable on a surface of constant y + , have K e R. For 
K < this corresponds to unstable modes with p + purely imaginary and y complex. 

Let us first focus on the case A < and, for simplicity, m — 0. The classical trajectories 
and behavior of wave functions can be immediately understood from figure 6. There are 
three different kinematic regimes, determined by the light-cone energy p + . For < p 2 < 
(wR) 2 the string comes from y = — oo until it reaches a turning point yo > 0, where it 
bounces back. The wave function oscillates up to the turning point and becomes evanescent 
thereafter. In the low energy limit p + = 0, the slope of potential remains positive and the 
wave function oscillates for y < y . This oscillatory behavior is analogous to a tachyon 
in Minkowski space, whose wave function also has spatial oscillations for low energies. 
We shall come back to this important issue later because these states will be crucial in 
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Figure 6: Effective potential V(y) and conserved energy —A for string motion along direction 
y. Shown is the case of A < and m = 0, for p 2 < (wR) 2 (left) and p\ > (wR) 2 (right). For 
p 2 = (wR) 2 the potential becomes flat. In contrast with the particle case, when p + = the 
potential is tilted. 



protecting chronology. As one increases p + , the slope decreases and vanishes for p 2 = 
(wR) 2 . At this point the string moves at constant velocity, reflecting a cancellation between 
the string tension and the inertial force due the accelerating y-frame. The corresponding 
wave function is a plane wave. Finally, when p 2 > (wR) 2 the string comes from y = +oo 
and turns around at yo < 0. Its wave function becomes evanescent beyond that point. 
For large light-cone energies, the potential becomes infinitely steep, as it was the case for 
particles. 

In figure 7 we plot the dispersion relation for states with A < and m = 0. As shown 
in appendix A. 2, to determine the evolution of generic perturbations denned on a surface 
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Figure 7: Dispersion relation for string states with A < and m = 0. 
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Figure 8: Effective potential V(y) and conserved energy —A for string motion along direction 
y. Shown is the case of A > and m = 0, for p 2 < (wR) 2 (left) and p 2 > (wR) 2 (right). For 
p 2 = (wR) 2 there is no classical trajectory nor wave function. In contrast with the particle case, 
when p + = the effective potential is tilted. 



of constant y + , one must include all real values of yo in the on-shell basis. Hence, we must 
consider states with imaginary p + , which grow in time. The turning point yo for these 
unstable states satisfies < y < y c ( w ) with 

vM = 2E^RJ 2 ■ 

One then expects that the condensation of these modes will drastically change spacetime 
in the region y < y , where the wave function oscillates (see figure 6, left). 

We move to the A > case, and again set m — 0. Trajectories and wave functions can 
be understood looking at figure 8. For < p\ < (wR) 2 the string comes from negative y 
and turns back at yo < 0. The low energy limit p + = is similar to the above case but 
with y Q < 0. As the potential tilts and one reaches p\ = (wR) 2 , the turning point y is 
pushed to — oo and the wave function disappears. For p 2 > (wR) 2 the string comes in 
from positive y and bounces back at yo > 0. The high energy limit is the same as before. 
In figure 9 we plot the dispersion relation for A > and m — 0, both for p + real and 
imaginary. The unstable modes will grow considerably for y < yo < 0. 

For a fixed value of A and generic quantum number m not all real values of yo give 
rise to normalizable wave functions on a surface of constant y + . The condition for having 
on-shell stable states with real yo can be determined from the quadratic equation (11) 

2Ey U + ^-y + w 2 R 2 {yo)\ > , (14) 
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Figure 9: Dispersion relation for string states with A > and m = 0. 

where R 2 (yo) = 2Ey R 2 is again the proper radius squared of the orbifold compact direc- 
tion, computed at the turning point yo- On-shell unstable states, which are normalizable 
on a surface of constant y + , occur for p + purely imaginary and yo complex fixed by the 
dispersion relation (11). 

3.2 Canonical quantization 

We now quantize the string using ^-coordinates and in light-cone gauge. The worldsheet 
theory with the Lagrangian (7) is not a free field theory. However, the worldsheet field Y~ 
obeys the free wave equation. Thus the light-cone gauge can be implemented by demanding 
that Y~ does not have oscillatory terms. In this gauge, the solutions of the wave equation 
for the worldsheet fields have the following expansion 

Y~{t,&) = 2Rwa + y~(T) 

Y(r, a) = y(r) + -±= £ - (a n e 2m ^) + a n e 2m ^) 

Y+(r,a) = y + {r) + ^Y,-{^)e^^ + aUr)e^^). 

As usual in light-cone gauge, the coefficients a+(r) and a+(r) are not independent degrees 
of freedom because they are fixed by the Virasoro constraints L n = L n = for n ^ 0. 
From the string Lagrangian (7) we obtain the conjugate momenta 

1 / • , \ 1 • 1 • 1 

p - = -7T- ( Y ~ 2EYY-) , P = —Y , P + = -— Y~ = -— y~ , 
2tt V J ' 2tt ' + 2tt 27r y ' 

which can be integrated along a to yield the string center of mass momentum derived 
in (9). The worldsheet energy-momentum tensor T±± = d±Y fJ, d±Y fl enables us to compute 
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the Virasoro zero-modes 

Lo = h Jo T - da = i v ^ + l Ey {wR)2 + ir + N ' 

where iV and N are the usual light-cone oscillator contributions. The Virasoro zero-modes 
sum and difference are 

L + L = ^p^ + % (u>i?) 2 + N + N , 
L - L = mw + N - N . 

Given any N and N, studying the classical dynamics of the string center of mass reduces 
to the effective Hamiltonian problem described in subsection 3.2. The corresponding con- 
served quantity A is fixed by setting L — L — 0. 

Canonical quantization of the string is done by imposing equal-time commutation re- 
lations. This gives the non-vanishing commutators for the string center of mass 

[y _ (r),p_(r)] = % , [y(r),p(r)} = i , [y + (r) , p+(r)] = % 

and for the oscillators 

Since the oscillators decouple from the zero modes the normal ordering for the Virasoro 
zero-modes is the usual one. For the bosonic string one has Lq = Lq = 1 and physical 
states satisfy the on-shell operator relation 

M 2 (y) = -p^ = 2Ey (wR) 2 + 2 (iV + N - 2) , 

and the level matching condition 

N — N = mw . 

The quantum numbers for the string center of mass wave functions <t> p+ ^ yo>m will then obey 
the on-shell relation 

A = 2p + ^ + 2Ey (pi - (wR) 2 ) = 2 (n + n - 2) . 

We conclude that for the bosonic string both positive and negative values of A are allowed. 
The corresponding behavior of the string center of mass wave functions was studied in the 
previous subsection, where we saw that states with negative values of A grow significantly 
in the bad region, as well as a bit inside the good region. 
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3.3 Partition function 



i 



We now compute the bosonic string partition function and find divergences associated with 
the condensation of physical states. We shall see these are all infrared divergences. In the 
canonical formalism one takes the trace over the Hilbert space of states 

IT T 2 

where q 2mT and r = T\ + ir 2 is the modular parameter of the Euclidean torus, restricted 
to the fundamental domain T . Since we are computing the trace in light-cone gauge, one 
integrates over all the string zero-modes but keeps only the physical transverse oscillators. 
Using the explicit form of the Virasoro zero-modes it follows that 

Z = [ — Tr exp hiiir^N — N — mw) + 2ixt 2 (-\paP^ - Ey (wR) 2 - N - N + 2) . 
jt r 2 1 \ 2 J _ 

In the expression for the Virasoro zero-modes the string center of mass momentum and the 
number operators commute. Therefore one can perform the trace over the oscillators to 
obtain the same result as for Minkowski space. Integrating over momenta in the spectator 
directions and choosing the basis \p + , yo, m) to perform the trace over the remaining zero- 
modes we arrive at the expression 

Z = L + V 23 f —Z (t) I -p-dyop(yo,p + ) exp ( - 2niT 1 mw + ttt 2 \ ) , (15) 

JT T 2 m w J Z7T 

where Z = (27r v /T^) _23 |?](r)| _48 and we recall that on the basis \p + ,y ,m) the operator H 
defined in (10) has eigenvalue —A given by (11). 

3.3.1 Infrared divergences 

We want to understand the divergences in the partition function coming from the infrared 
region r 2 — > +00. In particular, we are interested in investigating the contribution of each 
string state to such divergence. For this purpose we cut the ultraviolet region of modular 
integration by setting r 2 > 1, in order to integrate over T\ to enforce the level matching 
condition. Recalling the expansion of the Dedekind eta function r\{j) we have 

00 

Z, = (2n^- 2 y 23 (qq)- 1 ]T d n d h q n f , 

n,n=0 

where d n and dn represent the degeneracy of left and right moving states. As for the 
particle case we define the analytic continuation of the p + integral as p + — > ip + . The 
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infrared behavior of the partition function can then be studied from 



Z ~ % 




) 



exp -2tyt 2 ( -ip+ — + Ey (p+ + (wR) 2 ^ + n + n - 2 



where the sum is restricted by the level matching condition for physical states. The infrared 
divergences come from states in the Hilbert space such that the real part of the argument in 
the above exponential is positive. Hence, given (m, w, n, n) and negative y , the integral is 
infrared divergent for p + large enough. Let us then focus on the region of integration with y 
positive and look for further infrared divergences. Using the saddle point approximation 
we obtain the large r 2 behavior of the p + integral. Then one is left with a r 2 integral 
dominated by an exponential with argument 



As for the particle case, we see that the condition for the integral to be infrared convergent 
is equivalent to the condition for the existence of on-shell stable states with yo positive, as 
can be seen from (14). For A > 0, this condition holds for all on-shell states that oscillate 
in the good region, which have yo > 0. String states with A < 0, do not always satisfy the 
above condition and therefore give rise to infrared divergences. In particular, the state with 
m — 0, w — 1 and n = n = renders the partition function divergent for y < 2/(ER 2 ). 
Excluding the usual closed string tachyon, this is the state that condenses the furthest into 
the y > region. 

3.3.2 Ultraviolet finiteness 

Up to now we have cut the ultraviolet region of modular integration to understand the 
infrared origin of the divergences. In string theory one expects no further divergences 
arising from the ultraviolet. Let us analyze the partition function without cutting off the 
ultraviolet region. The divergences we will find correspond to those found in the previous 
analysis. 

Consider the expression (15) for the partition function. Performing the Wick rotation, 
the sum in Kaluza-Klein momentum can be Poisson re-summed using 
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such that the p + integral becomes trivial. As explained in appendix A.l, we can then 
use the integral representation of the measure to interchange the integral over y with the 
integral over the string center of mass y, arriving at the final expression 

where T = wt — w'. In appendix B we derive independently this result using the path 
integral formalism. It is easy to check that the integrand is modular invariant. Unlike the 
particle case, for fixed (w,w') and y — > there are no ultraviolet divergences because the 
r 2 integration is restricted to the fundamental domain T . 

3.3.3 Hagedorn behavior 

The term in the partition function with w = w' = gives the usual Minkowski space 
contribution. The remaining terms in the sum are the contribution from the winding 
states and Kaluza-Klein states in their dual guise. From now on we consider only these 
terms since they are the ones which contain the orbifold information. Using the usual trick 
of replacing the fundamental region T by the strip T = |r : — | < T\ < |,t 2 > Oj, while 
simultaneously setting w = and w' > in the sum [50], we can write 

If one sets R(y) = R this expression would be that of the partition function for bosonic 
strings in Minkowski space with one compact direction of radius R, which for a critical 
value of the moduli R has a divergence associated to a Hagedorn phase transition [3]. 
Similarly, in the O-plane orbifold we expect to obtain a divergence for a critical value of 
the string center of mass coordinate y. To see this expand the Dedekind eta function rj(r) 
in a series and perform the T\ integration to obtain 

z = i L + (2nR)V 23 ^ j2 fj f°°dr, 



(27T) 



2(i 



E dn fdy r ^ exp (-An(n - l)r 2 - ^ w 2 R 2 (y) 



Let us analyze the r 2 integral for specific values of n. For n — 0, the integral diverges 
exponentially as r 2 — > oo. This is the usual infrared tachyonic divergence of the bosonic 
closed string. Note that the y dependence is irrelevant in what concerns this divergence. 
For n > 1 fixed, the r 2 integral diverges for y — > 0. It diverges in the region r 2 — > 0, so one 
might think that this is an ultraviolet divergence. However, this is not the case because the 
small r 2 region of the strip T does not result only from a modular map of the fundamental 



20 



domain ultraviolet region. We expect these divergences to be infrared divergences. An 
instructive analogy is to compare with the case of bosonic string compactified on a circle 
in the limit of vanishing radius. In this limit there are winding states that become massless 
and the divergence should be thought in the T-dual picture as a large volume effect. 

For positive y, the r 2 integral converges and the partition function can be explicitly 
written as a Hankel function of first kind 

Using the asymptotic expansions for the degeneracy of states and for the Hankel function 

d n ~ n" 27 / 4 e 4 ^ , H?\z) N ~°° ./A e <(™/2--/4) 

V 

we can see that the large n expansion of the integrand in the partition function is dominated 
by the exponential term 

We conclude that, for each winding number w, the sum in n diverges when 

, \ 2 
y < y c \w) = — — —7 . 

y y v ; E(wR) 2 

This behavior of the partition function for large n is of the same type of the Hagedorn 
divergence of the bosonic string compactified on a circle, a divergence which can be traced 
back to the appearance of unstable tachyonic modes in the spectrum [4] . 

4 Hagedorn transition - the chronology warden 

Having found unstable states and related infrared divergences in the string partition func- 
tion that resemble what happens in the Hagedorn phase transition, we now suggest how 
this can protect chronology. 

Let us start by recalling the behavior of a tachyon field in Minkowski space. The 
condensation of the field occurs via the growth of low momentum modes, as can be seen 
from the dispersion relation plotted in figure 10. In the case of the bosonic string one 
expects that the closed string tachyon condensation will change the spacetime structure. 
For the O-plane orbifold a similar phenomenon occurs. In particular, there are winding 
string states that play an important role because the associated field has unstable modes 
that propagate not only in the bad region, but also a bit inside the good region. The growth 
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of these modes will then induce a large backreaction in the region where they propagate, 
triggering a phase transition. It is therefore natural to conjecture that this phase transition 
protects chronology. 

For simplicity, let us focus on the bosonic string and we shall comment on the extension 
to the superstring afterwards. As usual, we neglect the closed string tachyon. Consider 
the state with m — 0, w — 1 and n = n = 0. From equation (13) we can compute the 
position where the mass M(y) of this field vanishes 

y = yc = ^- (i6) 

For y < y c the mass squared of the field becomes negative. Hence we expect a phase 
transition in this region. This expectation is confirmed by the existence of on-shell states 
with imaginary light-cone energy. These states grow significantly for y < y , with the 
quantum number yo satisfying < yo < y c , as shown in figure 7. The singularity developed 
by the partition function for y < y c confirms these results. 

The O-plane orbifold phase transition shares some properties with the Hagedorn phase 
transition. Consider the bosonic string compactified on circle of radius R. The mass of a 
string with momentum and winding quantum numbers (m, w) is 

77? 2 

M 2 = 2 (n + n - 2) + — + (wR) 2 . (17) 

States with quantum numbers (0, 1) and no oscillator excitations have M 2 = — 4 + R 2 , 
becoming massless for R = 2. This is precisely the proper radius of the O-plane orbifold 
compact direction when the state with the same quantum numbers becomes massless. In 
the usual Hagedorn phase transition the radius R is a moduli. As the radius decreases to 
R = 2 the first winding state w — 1 condenses and triggers the phase transition. On the 
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other hand, for the O-plane orbifold the compactification radius varies with the coordinate 
y. Since the w — 1 winding state is the one that penetrates the furthest into the good 
region, one expects that the condensation of this field will cause significant changes on the 
spacetime structure for y < y c . In other words, the string fields will develop a large tadpole 
beyond this critical point. 

4.1 Extension to the superstring 

It is straightforward to extend the previous results to the superstring, where one needs to 
specify the orbifold spin structure. For both spin structures, we need to recall that, in the 
NS-NS sector, the normal ordering for the Virasoro zero-modes gives the on-shell relation 
A = 2(N + N-1). 

Consider first the supersymmetry breaking case where spacetime fermions have anti- 
periodic boundary conditions [51]. In this case states with even winding number have the 
usual supersymmetric GSO projection, while states with odd winding number have reversed 
GSO projection. The physical picture is very similar to the bosonic string just considered, 
but one does not have the usual bosonic string tachyon. There are winding states that 
condense and whose wave function also oscillates inside the region y > 0. These states 
have A = —2, m = and odd winding number w. For w — ±1 the mass M(y) vanishes 
for y = y c = 1/(ER 2 ) and it becomes negative for y < y c . Correspondingly, on-shell states 
with < yo < y c have imaginary p + and grow significantly for y < y . The proper radius of 
the compactification circle at the critical value y = y c is equal to the superstring Hagedorn 
critical radius. Moreover, the partition function has infrared divergences for y < y c . 

The case of supersymmetric boundary conditions is qualitatively different. After im- 
posing the standard GSO projection only states with A > survive. Now there are no 
winding states that condense inside the good region. The contribution from each string 
field to the partition function will then diverge for yo < 0. However, there will be the usual 
supersymmetry cancellation between bosons and fermions. All unstable states condense 
in the bad region. This is consistent with the conjecture put forward in [5] to protect 
chronology in the presence of closed causal curves, which we now recall. 

Start with the O-plane orbifold in type IIA theory and ignore the y < region. Next 
uplift the geometry to M-theory and reduce back to type IIA along the orbifold direction 
y~ . After performing a T-duality transformation along the eight spectator directions, the 
result is the geometry of an 08-plane orientifold [52]. This is the reason why this orbifold 
was called the O-plane orbifold. The above duality chain suggests that one should simply 
excise the y < region of the orbifold. Moreover, the quantization of the 08-plane charge 
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was seen to be dual to a quantization of the string coupling constant g s in the orbifold side. 
Then it was conjectured that such coupling quantization is necessary to restore unitarity in 
the O-plane orbifold, lost at arbitrary values of the coupling due to the presence of closed 
causal curves. 

The results of the present paper fit-in nicely with the above conjecture. Our computa- 
tions are valid at zero string coupling, where we see that the string fields condense in the 
region y < 0. It is therefore conceivable that the end-point of the supersymmetric O-plane 
orbifold condensation results in the excision of the bad region and in a vacuum expectation 
value for the dilaton field consistent with dual 08-plane charge quantization. 

5 Concluding remarks 

We believe the results here presented for the O-plane orbifold with supersymmetry break- 
ing boundary conditions are general and provide a truly stringy mechanism to protect 
chronology. Suppose we have an arbitrary spacetime with a Cauchy surface. Suppose also 
that, as the geometry evolves, we are about to create a closed null curve in the future of 
such surface. Then, when there are closed spacelike curves with proper radius of order of 
the string length, winding string fields become massless and start condensing. There will 
be a Hagedorn phase transition, which protects chronology. In other words, after the phase 
transition the new vacuum will be chronological. 

Hawking's chronology protection mechanism is based on an ultraviolet effect. He ad- 
vocates that strong curvature effects due to a large one loop energy-momentum tensor 
for the matter fields will create a singularity or will prevent spacetime to curve such that 
non-causal curves do not form. In contrast, the string theory mechanism for chronology 
protection is an infrared effect, as it is usually the case for phase transitions. 

Other situations in string theory where the appearance of light winding states triggers a 
phase transition are well known. An interesting follow up of our proposal is to find the exact 
end-point of the condensation. Our conjecture is that this phase transition will remove 
the bad region from spacetime. Indeed, phase transitions associated with light winding 
string states, that produce such a drastic effect in spacetime and that lead to a change of 
topology, are expected to occur in the context of the AdS/Gauge-theory duality at finite 
temperature [53] , or in the context of circle compactifications with supersymmetry breaking 
boundary conditions [54]. The authors of the latter work argued that, after condensation 
of the winding string field, the worldsheet theory will flow to an infrared trivial fixed point 
with no target space geometric interpretation. We believe something similar will happen 
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that excludes the formation of non-causal curves. 

In order to check the generality of the mechanism presented herein, interesting labora- 
tories include other time-dependent orbifolds and other closed string backgrounds, which 
have already been studied in relation to closed causal curves in string theory. 
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A Wave functions in the O-plane orbifold 



In this appendix we prove several results used in the main text regarding wave functions. 
In the appendix A.l we construct a complete basis of states in the orbifold geometry and 
define a measure for these states. In appendix A. 2 we consider on-shell wave functions 
normalizable on a surface of constant y + . 

A.l Off-shell wave functions and Hilbert space measure 

To define a measure in the Hilbert space we shall regularize the wave functions by working 
in a box of volume V = 2nRLL + . More concretely, we impose the boundary conditions 

(y~,y, y + ) = <t> (y~ + ^R, y, y + ) 
(y~, y, y + ) = <t> (y~,y, y + + L+) 

<t> [y~, -7p y + ) =<p(y~, ^, y + ) = o . 

Furthermore, it is convenient to work with a basis of functions formed by eigenfunctions 
of the operator 

H = 2<9 + <9_ + 2Ey (dl + (wR) 2 ) - d 2 y . 
These wave functions can be written in the form 

= f( z ) e i( P +y + +p-y-) ? 

where p + = 2nr/L + and p- = m/R, for r and m integers. For fixed winding number w, 
we define z 3 = K(y — y) 3 , with the constant K depending on the light-cone energy as 

K = K( P+ ) = 2E [p% - (wR) 2 ) 

and with y related to the eigenvalue —A of the operator H by 

A = 2p + p_ + y K(p+) . 

The function f(z) satisfies the Airy differential equation 

(d 2 - z) f(z) = , 

subject to Dirichlet boundary conditions at the points 

z 3 = K ^y ±r]^j , 
26 



where r\ = sgn(K) and we conveniently chose z + > Z-. With these boundary conditions 
the Airy differential equation has solutions only for discrete values of y , which we denote 
as 2/0 = yo(s), with s a positive integer. The spectrum of the operator A is now discrete. 

To define the Hilbert space measure in the large volume limit we need to understand the 
behavior of the function yo(s) for large L. Since the function f(z) is a linear combination 
of the Airy functions Ai(z) and Bi(z), the boundary condition f(z+) = gives 

f(z) oc Bi(z+) Ai(z) - Ai(z+) Bi(z) . 

To quantize yo one imposes the other boundary condition 

Bi(z+) Ai(z_) = M(z+) Bi(z_) . 

The density of states around a fixed value of yo can be computed by observing that in the 
limit of large L one has z + >> 1. Then the above equation reduces to Ai(z_) = because 
of the exponential growth of the function Bi(z) for positive real values of z. Denoting the 
zeros of the Airy function Ai by a s , the function yo(s) is implicitly defined by 



a s 3 = K (y (s) 



L 

V 2 



Again in the limit of large L, one can use the asymptotic expansion for the zeros of the 
Airy function 

2/3 



to determine yo(s). Two consecutive zeros will then give the following difference between 
discrete values of yo 

3tt(4s - 1) 1 " 1/3 



Syo = yo(s + l) - yo(s) ~ -it 



K 



which can be written as a function of yo 



-IT 



1-1/2 



Since 5y — > for L — > +oo we can move from sums to integrals 

£ F Ms)) - [ 2 dyop(yo)F(y ), 



with measure p = E s ^ (yo - 2/o(s)) - Vl^ol given by 



p(yo,p+) ^ 



7T 



1/2 
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Notice that there may be other states with y outside the above region of integration, but 
these are negligible for large L. 

The normalized basis of functions in the large volume limit is 



p+,yo,m 



(y) = 



11/3 



Ai(z) e 



>27rRL + p(y ,p + ) 
satisfying the orthogonality and completeness conditions 

<%o - y' ) 8 m>m > , 
dyop(yo,p+) 4>* P+ , m , m {y') <t> P +,m, yo {y) =$(y-y') ■ 



2tt 1 
(P+, yo, m\p' + , y' , m') = —S(p + - p' + ) — , 

^+ p{yo,p+) 



^ 2vr 



It is important to have an independent check of the above measure. In the reminder of 
this appendix we provide that check. We shall start from the canonical partition function 
for a particle and derive the result obtained independently in the next appendix using the 
path integral formalism. This involves a Poisson summation of the Kaluza-Klein quantum 
number m and a rewriting of the measure p(yo,p+) as an integral over a variable y. From 
the path integral computation one sees that this integration variable should be interpreted 
as the particle average position along the y-direction. A similar calculation can be done 
for the string partition function. 

We start with the canonical particle partition function written as 



z = *r i ? ^ si dp+ /4 dy ° p{y °' w+) exp 

and write the measure p(yo> W+) m the integral form 



- [M 2 -2ip + - + 2Ep 2 + y 



p{yo, ip+) = 



\K(tp 



+, 



I V2 rk 



2tt 



Q[y(y-yo)] 
>Jv(y - yo) 



where is the usual Heaviside step function. The integral over y reduces to a half 
Gaussian integral in the L — > oo limit 



dyo 



Q[v(y 



yo 



f2 +1 )y du 

\fn{y - yo) " '° ^ 

Using this result the partition function becomes 



e 2 



(Ky-\K\u) 



dy 



exp 



2tt 
l\K\ 



L(M 2 -2tp™+2Eply 



1/2 



e 2 



Ky 



R 
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Then, performing the Gaussian integral over the light-cone energy we have 



T ff , f°° dl 1 R ^ 
Z = tL + J dyj 



■ - M H 

2 1 R 2 (y)J 



21 2-kI R(y) 

where R 2 (y) = 2EyR 2 as defined in the main text. The Poisson summation formula 



e 



can be used to obtain the expression for the partition function derived from the path 
integral formalism 



r0 °dl 1 (I n/r2 (2nwR(y)f 



/* poo di \ ll 



/o 21 (2tt/) 3 /2 ^ ^ 2 2/ 
A. 2 On-shell wave functions 

We shall consider the light-cone evolution of a generic field configuration defined on a 
surface of constant y + . Then, given the field and its first derivative on this surface, the 
evolution is determined by the equations of motion. We shall address this problem by 
considering a basis of normalizable functions on the surface of constant y + with simple 
on-shell evolution. 

The basis can be constructed using the orbifold modes <f>p + , yo ,m(y), restricted to a con- 
stant y + surface and with the quantum numbers satisfying the on-shell condition 

X = 2p + p^ + 2Ey (pl-(wR) 2 ) , 

for A fixed. We shall call these functions 4> yo ,m{y~ ■> y)- Discarding an irrelevant normaliza- 
tion constant, we define these functions by 

KAViV) =K^M(K 1 / 3 (y -y)) . 

Their on-shell evolution to the orbifold spacetime is determined by the on-shell value of 
the light-cone energy p + . From the asymptotic behavior of the Airy functions, it follows 
that 

K = 2E (pi - (wR) 2 ) 

must be real, otherwise the function Ai would grow exponentially for y — > +oo or for 
y — > — oo. Notice that, one can not include the other Airy function Bi, because it always 
grows exponentially at least in one direction. We conclude that the functions <f> yo ,m(y~,y) 
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Figure 11: Contours in the /%o complex plane for \a\ < 1. The lines outside the real axis belong 
to an ellipse with foci at j± = — 1 ± y/l — a 2 . For \a\ = 1 the ellipse becomes a circle and the 
contours start, instead, at the center of the circle. For a = the ellipse degenerates to the real 
axis. For \a\ > 1 the foci of the ellipse move to the complex plane, the left contour starts and 
finishes at (3yo — > +oo and the right contour at (3yo — ► — oo. 



are normalizable provided both y and p + are real, or y is complex and p + is purely 
imaginary. 

The set of solutions to the on-shell condition with K real can be parameterized as 2 



(3X 
~K 



u 



/3yo = u ± a \Ju 2 + 2u , 



wR — — 

p+ = =F vr + 2u , 

u 



where «6l and we denned a = 2mw/\ and (3 = E(2wR) 2 /\. In this parameterization 
the unstable modes appear for — 2 < u < 0. For each value of u there are two functions 
<t>u m (y~ , y) , corresponding to the choice of signs in the above parameterization. In the 
general case m/0, these functions are different because they have different values of yo- 
In figure 11, both contours in the yo complex plane for on-shell states are shown. In the 
simpler case m — 0, both contours coalesce to the real line and, for each value of u, the 
two functions <f>u,m(v~ ■> v) are e Q ua l with symmetric values of the light-cone energy. 

Using the above parameterization, the initial value problem reduces to finding the 
functions F ± (m, m) such that 



F (y,y ) = / du [ F+ ( u ^ m )Km(y^y ) + F («, m ) K, m (y, y~ 

F '(y,y~) = -iwRj2 / du u~Wu 2 + 2u \P + (u, m) (p+ m (y, y~) - F~ (u, m) (p~ m (y, y 



2 The case of w = can be included by setting u — w 2 t and then taking the limit w — ► with fixed t. 
In this limit, the ellipse in figure 11, associated to complex values of j/Oj degenerates to a parabola. 
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for any given function F(y, y ) and its first derivative F'(y, y ). We start by decomposing 
the function and its derivative in Fourier modes 



F(y,y~)=J2 [dq e^ +i ^~ G m (q) , 
F'(y,y-) = J2 fdq e^ +i ^~ G' m {q) . 



Then, by linearity, we can break the problem in two cases: F generic and F' = 0; F' 
generic and F — 0. Moreover, we only need to consider each Fourier mode independently 
Firstly, we consider an initial configuration with zero derivative 

F(y,y-)=e^ + ^y- , F'(y,y') = 0. 

Using the integral representation of the Airy function 3 



<fcjm(y> y ) = P e l * y J 7^ ex P is (Py-uT a ^ u<2 + 2u ) - is 



3A 



we find that the corresponding coefficients F^{u,m) in the expansion have to satisfy 



2tt 



5(q-f3s) = 

= J du u-Wu 2 + 2u e" m ( s+£ ^) 
Secondly, we choose the initial condition 
F(y,y-)=0, 



—isa\/u 2 +2u 



F + (u,m) + e 



isa\J 'u 2 +2u 



F q (u,m) 



e- isaVu ' +2u F a + (u,m) -e 



isa\/ u 2 +2u 



F (u,m) 



F'(y,y-) = e^ + ^y- , 



and find the following equations for the corresponding F /=t (-u, m 



-I 



du e 



e- isaV ^ T2 ^F' + (u,m) + e 



isa ^+^ F'-(u,m) 



— u~ 1 \ / u 2 + 2u e 

Z7T 



X 



X 



-isoV u 2 +2u pi + ^ u m ^ _ e isaVu 2 +2u pi 



F' (u,m) 



3 This integral representation is valid for complex values of yo(u), provided one deforms the real axis 
integration contour as s — > s + iesgn(u). 
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Finally, the expansion for the generic function F with derivative F' has coefficients 

F ± («,m) = Jdq[G m (q)F q ± (u,m) + G' m (q)F' q ± (u,m) . 

Unfortunately, we do not know if the basis of functions is complete because we were not 
able to compute the functions F^(u,m) and F'^(u,m) in the general case. However, in 
the case of functions independent of y~ , so that we only need to consider wave functions 
with m = 0, we showed that the basis is complete for A > and that it spans very generic 
functions for A < 0, as explained below. The latter case is the most relevant case for 
the results of this paper, since it includes the winding states that were argued to protect 
chronology in the O-plane orbifold. 

For m = the on-shell relation becomes A = yoK, and therefore reality of K implies 
that yo is also real. In the u parameterization this corresponds to setting a = and the 
above formulas simplify considerably. In particular, defining 

F, + (n,0)=F-(u,0) = ^F,(u) , 
F'+(u,0) = -F'-(u.O) = ——*L==F,(u) , 
one is left to find F q such that 

f du -iu( s+4f- I ^ / \ 
6(q-0s) = J — e V ) Fq {u) . 

For A > this equation has solution 




showing the completeness of the basis for A > (and A = by adding a mass regulator). 
When A < 0, the function s + /5 2 s 3 /3A is not injective in the interval \(3s\ < 2V— A . 
Therefore, the above explicit solution for F q is only valid for \q\ > 2^—\ . For \q\ < 2^/— A , 
the integral equation for F q does not have a solution, showing that the basis is not complete. 
The subspace of functions generated by this basis is the set of functions whose Fourier 
transform Go(q) is a function of the particular combination q + q 3 /3\. Hence, the Fourier 
coefficients with \q\ > \J — 3A can be fixed arbitrarily, but also determine the remaining low 
momenta Fourier modes. We conclude that we can construct well localized perturbations, 
with some constraint on their long range decay properties (this can be deduced from the 
constraints on the pole structure of the Fourier transform). An interesting question is if 
one can complete the basis or, on the other hand, if the above restrictions are necessary 
for consistent evolution in the presence of closed causal curves. 
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B Path integral formalism 



In this appendix we compute the particle and string partition functions using the path 
integral formalism. The expressions we obtain can also be derived from the canonical 
formalism using the method described in the previous appendix. They serve as a check of 
the measure in the Hilbert space and provide a physical interpretation for the integration 
variable y. 

Before we start the computations it is convenient to write the orbifold identifications 
more compactly as 

Q{R, A)x = e~ 27rAJ x + 2ttR f (ttA) , 
where we defined the column vector T(a) and the nilpotent matrix J as 



T(a) = 



( 1 \ 

—a 
2 2 

v r ) 



( \ 
1 

v o i o ; 



A straightforward computation shows that the action of the orbifold group on a vector x 
satisfies 

x n = Q n (R, A)x = Q(nR, nA) x . 
B.l Particle partition function 

The one-loop vacuum energy for a field of mass M is given by the particle partition function 



00 dl 
21 



VX exp 



_! M 2 -i f X 2 dr 
2 21 Jo 



where the functional integral is done over periodic trajectories. Working in covering space, 
these trajectories are periodic up to the orbifold identification 

X( T + l) = e- 27TwAJ X{r) + 2nwRf (nwA) , 

where w is the trajectory winding number. A general path can then be expanded in a 
Fourier series as 



X{t) = e 



-2-kwAtJ 



^X„e 2 ™ r + 2irwRrf (irwAr) , 



where X n = Xl n so that X is real. Replacing this expansion in the action the contributions 
from zero-modes and from quantum fluctuations decouple, so the partition function has 
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the form 



The contribution from quantum fluctuations is given by 4 



■'qu 



J (jj 2 3 ^X n rfX n ^ exp (- f; X n tM n X n ^) 



fr (2vr) 3 
iii detM„ ' 



where the matrix (M n ) MJ/ is 



M B = 



(27TO) J 





2wA . 






1 


-1 




n 











1 


1 


n 






-1 









o J 

The determinant is therefore equal to that in the Minkowski space computation and one 
obtains after zeta function regularization Z qu = i(2-7r/)~ 3 / 2 . 
Next we consider the contribution from the zero-modes 



dX exp 



(2tiwR) 2 E 
I 



(*o + § {X ) 



where the integration is over the orbifold fundamental domain. Moving to the ^-coordinates 
this reduces to 



-cl 



J dy dydy 



exp 



(2nwR(y)Y 
21 



Integrating over the y + and y directions, we obtain the final answer for the particle 
partition function 



/oo r 
dy / 
-oo JO 



-°4_4_ exp f-i M --< 2 ^)'" 

21 (2vr/) 3 / 2 \ 2 



21 



B.2 String partition function 

We now compute the bosonic string partition function 

Z = f OT2>ye" s ^ 7 l 



4 The factor of 2 3 in the measure appears because the coefficients associated to real normalized modes 
are y/2Re(X n ) and y/21m(X n ). 
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where the Euclidean worldsheet action is 



S = ~L I d2,J V71 ab daX ■ d b X 



The dot denotes the inner product A-B = A T GB, with G the target space metric. Using the 
symmetries of the action we fix the metric on the torus 7^ to be given by ds 2 = \da\+rda 2 | 2 , 
with modular parameter r = T\ + ir 2 . This gauge fixing introduces ghost contributions 
which will be added in the end. With this choice of metric, and following the notation 
introduced above, the orbifold identifications defining the (w±, w 2 ) -twisted sector are 

X{a x + 7T, a 2 ) = e- 2 ™ lA ^X(a 1 , cr 2 ) + 2nRw 1 f(nw 1 A) , 
X(a u a 2 + n) = e~ 27rW2AJ X(a ± , a 2 ) + 2nRw 2 f (nw 2 A) . 

Functions that obey these conditions can be written as a Fourier series 

X = e~ 2AvJ Xn U n 2 e 2i{niai+n2a2) + 2Rvf(Av) , 

n\,n 2 

— » — » — » — * 

where v = w\ai + w 2 a 2 and reality of X requires X_„ li _ ri2 = X* t n2 . The X 0i o are zero- 

— * 

modes while the other X ni ^ n2 describe quantum fluctuations. 

The gauge fixed partition function without the ghost contribution has the explicit form 

Z = J^d 2 r J VX exp J d 2 o ^ (|r| 2 <9iX • d x X - 2r x d x X • d 2 X + d 2 X ■ d 2 X) J , 

where the r integration is over the fundamental domain of the torus T . To evaluate 
the action of the above Fourier expansion of the fields it is convenient to consider the 
zero-modes and quantum fluctuations separately. We include the term 2RvT(Av) in the 
zero-mode contribution to the partition function. Start with the quantum fluctuations, for 
which we follow the analysis of [11, 16]. In this case the action has contributions of the 
type 

d a X ■ d b X = Y { 2m a C~ 2AVJ ~ ^ a AJ e- 2AvJ ) X nun2 e nma 1+ n 2 a 2 ) 

ni ,n 2 y^Q,0 

• Y fen b e- 2AvJ - 2w b AJe~ 2AvJ ) X nu1l2 e 2 ^+"^) 

ni,n2^0,0 

The structure of the various terms is such that after integrating over ai and a 2 we are left 
only with terms which, for each mode (rii,n 2 ), can be written as 

Xn u n 2 Crsi^J 1 ) T G J S X ni ,n 2 ■ 

\r,s=0 
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The path integral requires us to do the integration over the Fourier modes. These are 
Gaussian integrals and for each mode can be obtained from the determinant of a matrix 
(M niin2 )^ of the form 



M nuna = £ C rs (J r ) T GJ s 

r,s=0 



I M M ,. M_+ \ 

M x _ M xx 
V M+_ J 



It is then easy to see that the terms that give a contribution to the determinant of M ni>n2 
have r = s = 0. These terms only appear for w\ = w 2 = 0, so quantum fluctuations 
are insensitive to the winding of the string. The same happened in the computation of 
the partition functions for the null-boost and null-brane orbifolds presented in [11, 16]. 
Including the contribution of ghosts, the result is then the same as in Minkowski space, 
namely 

Z qu (r) = - (2n^r"\r ] (r)r . 

This is just the partition function of three free bosons plus ghosts. 

Next we evaluate the zero-mode contribution to the partition function. When n\ = 

— » 

n-2 = 0, X simplifies to 

X = e~ 2AvJ x + 2Rvf(Av) , 

where we defined X 0i o = x. The string center of mass X does not coincide with the zero- 
mode x. However, changing to the ^-coordinates, one can check that the y component of 
the zero-mode coincides with the string center of mass component Y . This is an important 
point since, in the main text, we interpreted the integration variable y as the string center 
of mass. A straightforward computation shows that 

d a X ■ d b X = 8ER 2 w a w b (x +^(x-) 2 ^= 8ER 2 w a w b y , 
which yields an action 

„ 2irER 2 - 

S = y TT , T = w\T — u>2 • 

Adding the oscillators contribution from the transverse directions, the full partition func- 
tion is then 



i 



L + (2nR)V 23 f d 2 r 48 _ r°° ( ir 2 -\ 
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